Griffiths phases on complex networks 
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Quenched disorder is known to play a relevant role in dynamical processes and phase transitions. 
Its effects on the dynamics of complex networks have hardly been studied. Aimed at filling this gap, 
we analyze the Contact Process, i.e. the simplest propagation model, with quenched disorder on 
complex networks. We find Griffiths phases and other rare region effects, leading rather generically 
to anomalously slow (algebraic, logarithmic, ...) relaxation, on Erdos-Renyi networks. Similar effects 
are predicted to exist for other topologies with a finite percolation threshold. More surprisingly, we 
find that Griffiths phases can also emerge in the absence of quenched disorder, as a consequence of 
topological heterogeneity in networks with finite topological dimension. These results have a broad 
spectrum of implications for propagation phenomena and other dynamical processes on networks. 
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Networks have become a paradigm in the study of com- 
plex systems After initial efforts devoted to uncover 
their non-trivial topological features, the focus shifted 
to dynamical processes occurring on them Q . Models of 
epidemics Q , such as the susceptible -infected -recovered 
(SIR), the susceptible -infected -susceptible (SIS), or the 
contact process (CP), have played a prominent role 
The concept of "epidemics" covers a broad variety of 
processes including real epidemics, computer viruses, ru- 
mor spreading, or signal propagation in neural nets. A 
remarkable finding is the absence of a finite infection 
threshold for the SIS process in heterogeneous scale-free 
networks; i.e. in contrast with regular lattices, where 
a critical point separates an active from an absorbing 
phase [H, these networks host endemic states even for 
arbitrarily small infection rates Q ; i.e. topological disor- 
der crucially affects dynamical processes on networks. 

Quenched disorder is well known to induce novel be- 
havior (such as, universality changes) in phase tran- 
sitions both in equilibrium and away from it. Un- 
der some circumstances, it may also generate phases 
unheard-of in pure systems. Consider, for illustration, 
the quenched contact process (QCP) (4|, characterized 
by a site-dependent quenched infection rate, on a lattice. 
Even if the system is globally in its absorbing phase, dis- 
order fluctuations can generate rare active regions with 
over-average infection rate. In these, activity lingers for 
extremely long periods; however, being finite, they in- 
eluctably end up falling into the absorbing state. The 
convolution of exponentially rare regions with exponen- 
tially large surviving times gives rise (see below) to a 
region in the absorbing phase, characterized by a generic 
algebraic decay of activity, i.e. a Griffiths phase (GP) 



This is just an example of a generic phenomenon, thor- 
oughly studied in the disordered systems literature in 
classical, quantum, and non-equilibrium systems oc- 



curring whenever exponentially distributed cluster sizes, 
P(s) ~ exp(— cs), have exponentially long activity times, 
t(s) ~ exp(6s); this leads to a GP with algebraic decay 
~ r _c / h with a continuously varying exponent (see 
for analytical approaches, further details and appli- 
cations, and Q for a recent review). Observe that, on 
the opposite to usual critical points, in GPs one observes 
scale-invar iance generically, i.e. without the need of pa- 
rameter fine-tuning. Similarly, other size distributions 
(as power-laws) do lead to different functional forms of 
slow relaxation. Let us stress that many disorder effects 
appear rather similarly in equilibrium and away from it; 
for instance, the (strong disorder) critical point of the 
(non-equilibrium) QCP is in the same universality class 
as the (equilibrium) random transverse-field Ising model 
0; i.e. rare-region effects transcend the frontier between 
equilibrium and non-equilibrium. 

From this broad perspective, it is surprising that very 
little attention has been paid so far to the effect of 
quenched disorder on dynamical processes on complex 
networks. Heterogeneity in the intrinsic properties of 
nodes is a very natural (not to say, unavoidable) fea- 
ture of real networks: node-dependent rates appear in 
all the examples above, owing to the specificity of the 
individual immune response, presence of antivirus soft- 
ware, and so forth. Our aim here is twofold: (i) we 
tackle the study of quenched node disorder in dynam- 
ical processes on complex networks, and look for rare- 
region effects in the simplest possible epidemic model, 
i.e. the QCP on Erdos-Renyi (ER) random networks; 
(ii) we explore whether network topological disorder on 
its own can induce rare-region effects. We report on the 
existence of Griffiths effects, including various non-trivial 
regimes with generic slow decay of activity, for different 
networks with quenched node disorder and/or topologi- 
cal disorder. Our main conclusion is that disorder, either 
quenched or topological, may induce under the conditions 
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studied here, slow relaxation on network dynamics; this 
is expected to go beyond the considered examples, and 
to apply to different models, dynamics and topologies. 

In the pure Contact Process, infected/active nodes heal 
at rate p and infect a randomly chosen nearest-neighbor 
(provided it was un-infected) at rate A. A simple (one- 
site) mean-field treatment of the average activity density, 
p, predicts an absorbing phase transition at A* = p, i.e. 
where the infection and healing rate compensate each 
other [8]. This prediction for the threshold value is not 
correct for finitely connected networks since in that case 
activity appears in localized regions, reducing the effec- 
tive infection rate (i.e. the prob. of having an already 
infected neighbor is larger than in the perfect mean-field 
mixing), i.e. A c > A*. A lengthy but standard [3] two-site 
or "pair" approximation, leads to A^ = p(k)/({k) — 1), 
where (k) is the average network degree |9(. Observe 
that the correction factor with respect to Aj converges to 

I when (k) — > oo (i.e. where mean-field holds) and di- 
verges at the percolation threshold (k) = 1 below which 
the network is fragmented @, [l(| , it cannot sustain ac- 
tivity and the phase transition disappears. 

We consider a QCP [4| with p = 1 and a quenched 
disordered infection rate: a fraction 1 — q of the nodes 
(type-I) take value A and the remaining fraction q (type- 

II nodes) take a reduced value Ar, with < r < 1. 
Obviously, for q — and q = 1 the model is pure, 
X c (q = 1) = X c (q = 0)/r, while for < q < 1 one 
expects A c to interpolate between these limits. For sim- 
plicity, we fix r = from now on. The average activ- 
ity density is p = (1 — q)p\ + qp2 (subscripts stand- 
ing for node type) and the one-site mean-field equa- 
tions are p\{t) = -p t + A(l - pi - /02)((l - q)pi) for 
i = 1,2, whose stationary solution has a critical point at 
Xl(q) = 1/(1 — q). As above, the one-site result needs 
to be corrected by a factor (k)/((k) — 1) to account for 
activity clustering, leading to: 
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(fc)-l l-q 



(1) 



Relevant for what follows is that type-I nodes experience 
a percolation transition where the type I-to-type I aver- 
age degree is 1 @, i.e. at q perc = 1 — (k)^ 1 . For 
Q > Qperc activity cannot be sustained: type-I clusters 
are finite and type-II ones do not propagate activity. 

We now present the results of a numerical investiga- 
tion of the QCP on ER networks with (k) = 3 (im- 
plying q perc = 2/3), and sizes up to N — 10 7 . Simu- 
lations are performed in a standard way: all sites are 
declared active initially and the dynamics proceeds as 
follows [3]. A site, i, is randomly selected and it ei- 
ther heals (with prob. 1/(1 + A;)) or infects a ran- 
domly selected neighbor provided it was empty (with 
prob. Ai/(1 + A;), where A; = A or depending on 
site type). We monitor the activity decay averaged over 
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FIG. 1: Average activity density p(t) vs time t for ER net- 
works with (k) = 3, r = 0, and N = 10 5 . As are ordered from 
top to bottom in all panels, (a) Upper panel: q — 0.6, and 
A = 5, 3.8, 3.6, 3.55, 3.5, 3.3; the dashed line is proportional to 
t -1 . (a) Inset: p vs ln(t) for q = 2/3; A = 10, 7, 5, 4.5, 4; the 
dashed line is proportional to ln(i) -1 / 2 . (b) Lower panel: 
q = 0.9, and A = 50,30,20,15,10,9,7,5,4.5,2.7. Straight 
lines lie in the Griffiths phase. 



many runs. Numerical results are synthesized in Fig. Q] 
and the inferred phase-diagram is summarized in Fig. [2] 
For q < 2/3 a critical line X c (q), which is very well fitted 
by Eq. (JT]) (implying that the pair approximation cap- 
tures the most relevant correlations), separates an ac- 
tive phase from an absorbing one (Fig. QJi). Instead for 
the fragmented case, q > 2/3 there is no active phase, 
as predicted above (Fig. [lb). The absorbing phase can 
be divided into various sub-phases. For q > 2/3 and 
X(q) > X c (q = q P erc) ~ 4.5 we observe (Fig.QJ)) a power- 
law decay with continuously varying exponents: 
p{t) - t-"^ x \ with 6 -> q -> q = 2/3 i.e. a CP. Right 
at q = 2/3 we find a logarithmic decay pit) ~ ln(t) -1 ' 2 
for A > X c {q perc ) w 4.5 (Fig. [Th, (inset)). For q > 2/3 
and X c (q = 0) < A < X c (q perc ), as well as for q < 2/3 
and X c (q = 0) < A < X c (q) (see Fig. [2]), we can fit a 
stretched exponential (not shown). Finally, below the 
threshold of the pure system X c (q = 0) ~ 1.5 the decay 
is purely exponential. 

Let us now rationalize these observations by using, as 
customarily done in disordered systems (see [3-lrj|. flljl). 
optimal fluctuation arguments. The following regimes 
can be predicted: (i) Griffiths phase: X > X c (q perc ) and 
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FIG. 2: Phase diagram for r — 0. See main text for details. 



q > qperc- The cluster size distribution for ER networks 



in the fragmented case is (12|: 
1 



P(s) 



2irp 



s -3/2 e -s(p-l-ln(p)) 



(2) 



where s is the cluster size and p is the average number of 
links per node, which in our case is p = (fc) 9= o(l — q) 
for type-I nodes. Within any given cluster of type- 
I nodes, let us define pi oc as the local average num- 
ber of links per node; and from it, a local value of 
q, qioc = 1 - Ptoc/{k) q =o- Obviously, for connected 
type-I clusters are above the percolation threshold, i.e. 
qioc < qperc, and hence, if A > A c (g perc ), they are active 
rare regions where activity survives until a coherent fluc- 
tuation extinguishes it. This occurs after a characteristic 
time r(s) which grows exponentially with cluster size, i.e. 
t ~ to exp(^4(A)s), where to and A(X) do not depend on 
s. Hence, the overall activity decays as: 



p(t)~ da S P(s)e X p[-V(*oe A(A)s )]. 



(3) 



Plugging Eq. @ into Eq. ^ and using a saddle point ap- 
proximation, one obtains p(t) ~ t~ e ( p ' X \ with 9(p,\) = 
— (p— 1 — ln(p))/yl(A); i.e. there is a generic power-law 
decay with continuously varying exponents, i.e. a 
GP. Note that 9 — > when p — > 1 , suggesting a crossover 
to logarithmic decay, (ii) Indeed, right at the percolation 
threshold, q = q pe rc, the cluster size distribution, Eq. ©, 
becomes a power-law leading, when plugged in Eq. ([3]), 
to a logarithmic decay, /?(t) ~ [ln(i/to)] _1 ^ 2 , for any A 
larger than A c (g perc ) (for which rare regions are active), 
(iii) For q < q perc there is a giant component of type- 
I nodes and, hence, standard mean-field like contact- 
process behavior is expected; i.e. power-law decay, 
p(t) ~ i -1 at criticality. Similar arguments lead to: (iv) 
In the absorbing region with X c (q = 0) < A < A c (g perc ), 
rare regions with qi oc < q perc exist but, in contrast to 
previous cases, they are subcritical (activity decays ex- 
ponentially fast in each of them); a stretched exponen- 
tial fits the decay. Active rare clusters can also appear; 



but they require qi oc to be smaller that the critical value 
of q for the given A, a stringent condition, leading (in all 
the region (iv) in Fig.©) to weak rare-region effects (as 
opposed to the previous strong effects) Q . Finally, region 
(v) with A < X c (q = 0) is free from rare-region effects: 
all clusters are subcritical with exponential decay. 

Summing up, optimal fluctuation arguments explain 
all numerical findings. Rare regions play a key role giving 
rise to generic slow decay of activity. Similar results hold 
for any value of r, i.e. when activity flows out of any site. 
Hence, strong rare-region effects are expected to occur 
in generic networks with a finite percolation threshold 
(e.g. structured scale-free networks [13[ or networks with 
communities) , while if the percolation threshold vanishes 
(e.g. Barabasi- Albert scale- free nets [2]) only weak rare- 
region effects are predicted. 

In what follows, we investigate whether GPs induced 
merely by topological disorder can exist. Let us first 
recall that Eq. ([TJ) provides an excellent approximation 
to the ((fc)-dependent) critical point. As a consequence, 
networks with heterogeneous degree may, in principle, 
exhibit topological rare regions: clusters with local de- 
gree, fc, above average would have a smaller percolation 
threshold and could be locally active even if the whole 
network is absorbing. For argument's sake, let us con- 
sider the CP on a network with bimodal degree distribu- 
tion, P(k) = p5(k - fci) + (1 - p)S(k - k 2 ) with kt > fc 2 ; 
a priori, one could expect rare active regions (with over- 
density of fci-nodes) to exist. However, numerically we 
find just conventional, non-disordered, exponential de- 
cay. Why is it so? In d-dimcnsional lattices disorder is 
known to be irrelevant for large d; the number of nearest 
neighbors is so high that the central limit theorem pre- 
cludes rare regions (i.e. deviations from normality) from 
existing. An extension of the concept of Euclidean di- 
mension to arbitrary graphs is the topological dimension, 
D, which measures how the number of nodes in a neigh- 
borhood grows as a function of the topological distance 
from an arbitrary origin: N(l) ~ l D [10(. In parallel with 
the result for lattices, we conjecture that GPs and similar 
rare-region effects do not typically exist for the CP in net- 
works with an infinite topological dimension (such as ER 
graphs above their percolation threshold or the discussed 
bimodal graphs), justifying our numerical findings. In- 
stead, for ER below the percolation threshold, which are 
characterized by a vanishing effective topological dimen- 
sion, D = (i.e. the number of nodes in any neighbor- 
hood/cluster converges to a constant for large values of I) 
GPs can exist. This suggests to explore CP on networks 
with finite D. As an example, we take the generalized 
random small-world networks defined in [141 ] as follows. 
Starting with a ring of L nodes, one defines the distance 
between nodes i and j as I = min(|i — j\,L — \i — all 
nearest neighbors are connected and any pair with I > 1 
is connected with probability P(l) = 1 — exp(— (3l~ a ). 
We have identified three different cases: i) For a < 2 



4 




FIG. 3: Density decay in generalized small-world networks 
with a = 2 and f3 — 0.2 for different values of A (from 
top to bottom: 2.81,2.795,2.782,2.77,2.75,2.73,2.71,2.70, 
2.69, 2.67, 2.65, 2.6). Straight lines lie in the Griffiths phase. 



the network diameter grows poly-logarithmically with N, 
hence, formally D = oo; indeed, as conjectured, no Grif- 
fiths phases are observed in our numerics (not shown). 
The critical behavior seems to correspond to that of 
the contact process with Levy- flight jumps |3|]. ii) For 
a > 2, (I) is finite, and correspondingly D = 1; long- 
edges do not alter the system dimensionality but they 
do introduce quenched disorder and a GP. We conjec- 
ture the critical behavior to coincide with that of the 
one-dimensional infinite randomness fixed point of the 
quenched contact process @, @] , as indeed verified numer- 
ically (not shown) ; topological and quenched disorder can 
lead to the same universal behavior, iii) In the marginal 
case a = 2, D has been conjectured to be /3-dependent, 
D = D(j3) 0, El]. We have checked that D is indeed 
an increasing function of /3 (13(0) = 1, D(0.2) « 1.21, 
13(0.5) fa 2.3). We expect the critical point A c (/3) to 
decrease when /3 grows (i.e. when (k) increases; see 
Eq-CQ)). Hence, for a given /3, the model must be in the 
active phase if A > A c (0) = 3.297848(22) (critical point 
in the one-dimensional lattice Q) and must be inactive 
if A < lim^oo A c (/3) = 1 (critical point for the fully con- 
nected graph). Therefore, A c (oo) < A c (/3) < A c (0) and 
the possible GP is bounded by this interval. Numerical 
simulations for /3 = 0.2 confirm (Fig. [3]) the existence of 
a GP with generic power-law decay for A in [2.65,2.81]. 
The width of the GP decreases with increasing /3 (i.e. 
the larger D the smaller the rare region effects); actually, 
preliminary results, suggests the existence of a finite up- 
per critical dimension above which the GP disappears |9j . 
We have found similar topological GPs for other small- 
world nets, and we predict them to emerge in many other 
networks, such as spatially embedded networks or fractal 
flfjj scale-free ones, with finite topological dimension. 

In summary, taking the contact process as an example we 
have shown that quenched disorder can induce GPs and 



other rare-region effects leading to generic slow relaxation 
for dynamical processes on ER networks. Similar effects 
are argued to appear on other topologies with a finite 
percolation threshold. Furthermore, heterogeneity in the 
topology may suffice to generate GPs on its own. This 
can occur only if the network topological dimension is 
not infinite allowing for rare regions to exist. Our results 
are expected to apply to dynamical processes other than 
the CP (even if each case should be carefully examined) . 
An inspiring application, illustrating Griffiths effects at 
play, is provided by a recent work reporting on algebraic 
"forgetting" times (reproducing experimental findings) 
in a simple model of memory 171 ] : this generic slow de- 
cay crucially depends on the neural network topological 
disorder. Social networks with heterogeneous communi- 
ties have also reported to exhibit generic slow decay and, 
what we interpret as, severe rare-region effects [181 ] . 
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